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The decay ipi'^S) KgKl is observed for the first time using ip{2S) data collected with the 
Beijing Spectrometer (BESII) at the Beijing Electron Positron Collider (BEPC); the branching ratio 
is determined to be B{ip{2S) K%Kl) = (5.24±0.47±0.48) x 10"^ Compared with J/V' ^ K%Kl, 
the tp{2S) branching ratio is enhanced relative to the prediction of the perturbative QCD "12%" 
rule. The result, together with the branching ratios of tp{2S) decays to other pseudoscalar meson 
pairs (tt^tt" and K^K~), is used to investigate the relative phase between the three-gluon and the 
one-photon annihilation amplitudes of ip{2S) decays. 



2 



PACS numbers: 13.25.Gv, 12.38.Qk, 14.40.Gx 



It has been determined that for many two-body exclu- 
sive J/V' decays, hke vector pseudoscalar (VP), vector 
vector (VV), and pseudoscalar pseudoscalar (PP) meson 
decays and nucleon anti-nucleon (NN) decays, the rela- 
tive phases between the three-gluon and the one-photon 
annihilation amplitudes are near 90° ^'or 
'0(25') decays, the available information about the phase 
is much more limited because there are fewer experimen- 
tal measurements. It has been argued that the relative 
phases in tMp^S) decays should be similar to those in J /ip 
decays 0,13, but the analysis of ipi'^S) to VP decays in 
Ref. 0] indicates this phase is likely to be around 180°. 
Another analysis of this mode though shows the relative 
phase observed in J/tp decays could also fit these de- 
cays 0, but it could not rule out the 180° possibihty due 
to the big uncertainties in the experimental data. There- 
fore it is important to measure phases in other '>i){2S) 
decay modes. 

In ip{2S) PP, the currently available measurements 
on 7r+7r~ and K~^K~ are from DASP 0, with huge er- 
rors. These do not provide enough information to extract 
the phase since there are three free parameters in the 
parametrization of the PP amplitudes 4, 5, 1(J|. The 
result for ip{2S) — > KgK^ is also needed to get the 
phase (lH. 

Furthermore, there is a longstanding "p7r puzzle" be- 
tween J/-!/; and ^{2S) decays in some decay modes; many 
'0(25') decay channels compared with the corresponding 
J/f/' decays are suppressed relative to the perturbative 
QCD predicted "12% rule" [iJi. It is of great interest to 
check this in more channels. 

In this letter, the first observation of "0(25) RgK^ 
is reported, and its branching ratio is used to determine 
the phase between the three-gluon and one-photon an- 
nihilation amplitudes and to test the "12% rule" be- 
tween J/0 and 0'(25) decays. The data used for the 
analysis are taken with the Beijing Spectrometer (BE- 
SII) detector at the Beijing Electron-Positron Collider 
(BEPC) storage ring at a center-of-mass energy corre- 
sponding to M^(2S)- The data sample contains a total 
of 14(1 ± 5%) X 10^ V'(25) decays, as determined from 
inclusive tp{2S) hadronic decays |l3 |. 

EES is a conventional solenoidal magnet detector that 
is described in detail in Ref . 11411 . BESII is the upgraded 
version of the BES detector [l5|. A 12-layer vertex cham- 
ber (VC) surrounding the beam pipe provides trigger in- 
formation. A forty-layer main drift chamber (MDC), lo- 
cated radially outside the VC, provides trajectory and 
energy loss (dE/dx) information for charged tracks over 
85% of the total solid angle. The momentum resolution is 
ap/p = 0.017v/l (p in GeV/c), and the dE/dx reso- 
lution for hadron tracks is ^ 8%. An array of 48 scintilla- 
tion counters surrounding the MDC measures the timc- 



of-flight (TOE) of charged tracks with a resolution of 
~ 200 ps for hadrons. Radially outside the TOE system 
is a 12 radiation length, lead-gas barrel shower counter 
(BSC). This measures the energies of electrons and pho- 
tons over ^ 80% of the total solid angle with an energy 
resolution of as/E = 22%/Ve (E in GeV). Outside of 
the solenoidal coil, which provides a 0.4 Tesla magnetic 
field over the tracking volume, is an iron flux return that 
is instrumented with three double layers of counters that 
identify muons with momentum greater than 0.5 GeV/c. 

A Monte Carlo simulation is used for the determina- 
tion of the detection efficiency. Eor the signal channel, 
■0(25) K^Kl, the angular distribution of the or 
is generated as sin^ 6*, where 9 is the polar angle in 
the laboratory system. The is allowed to decay and 
interact with the material in the detector, and for the 
K^, only — > tt+tt^ is generated. The detector re- 
sponse is simulated using a Geant3 based Monte Carlo 
program, SIMBES. Reasonable agreement between data 
and Monte Carlo simulation has been observed in vari- 
ous channels tested, including e+e~ — s- {•^)e^e~ , e+e~ — s- 
(7)/i+/i~, J/0 PP and ijj{2S) — > tt+tt^ J/0, J/-0 
l+l- {l^e,n). 

Candidate events are required to satisfy the following 
selection criteria: 

1 . The number of charged tracks is required to be two 
with net charge zero. 

2. Each track should satisfy | cos 6*1 < 0.8, where 9 is 
the polar angle in the MDC, and should have good 
helix fit so that the error matrix from track fitting 
is available for secondary vertex finding. 

3. £;*°* < 1.0 GeV, where is the total energy 
of the neutral clusters in the BSC which are not 
associated with the charged tracks. 

The two tracks are assumed to be tt"*" and 7r~. To 
find the intersection of the two tracks, an iterative, non- 
linear least squares technique is used. The intersection 
is taken as the vertex, and the momentum of the 
is calculated at this point. Eigure ^ shows a scat- 
ter plot of the TT^TT^ invariant mass versus the decay 
length in the transverse plane {L^y) for events that sat- 
isfy the above selection criteria and have tt+tt" momen- 
tum greater than 1.7 GeV/c. The cluster of events with 
mass consistent with the nominal Kg mass and with a 
long decay length indicates a clear Kg signal. The lack 
of events at L^y > 0.1 m is due to the trigger, which will 
be discussed later. 

Eits of the tt+tt^ invariant mass distributions (not 
shown) indicate good agreement between data and Monte 
Carlo simulation in both mass and the mass resolution. 
After requiring Lxy > 1 cm and 7r"'"7r~ mass within twice 
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FIG. 1: Scatter plot of tt n~ invariant mass versus the decay 
length in the transverse plane for events with ■k'^tv~ momen- 
tum greater than 1.7 GeV/c for a) data and b) Monte Carlo 
simulation. 



the mass resolution around the Kg nominal mass, 



the 



Kg momentum distribution, shown in Figure [3 is ob- 
tained. In the plot, there is a clear peak at around 
1.77 GeV/c with low background, as indicated by the Kg 
mass side band (three sigma away from the Kg nominal 
mass on both sides) events. The excess at lower momen- 
tum, which is not explained by the side band background, 
is due to the contribution of the background channels 
with Kg production. 

The main Kg production backgrounds near the sig- 
nal region are from il>{2S) —> K*°{892)K° + c.c. and 
i^{2S) 7Xci,Xci ^ K*° {892)1^ + c.c. where the K*° 
decays into K'^ and tt" and one of the K'^s becomes a 
Kg and the other one becomes a K^l. The background 
from -0(25") ^ 7Xcj(J = 0,2), Xcj ~^ K^K^ where one 
of the K^ decays into tt^tt~ and the other decays into 
■n^TT^ is removed almost entirely (more than 95%) by the 
cut, according to the Monte Carlo simulation. The 
decay ^{2S) J/^ + X with J/-0 decaying into K^X 
has a big branching ratio, but the K^ momentum is much 
lower. The light shaded histogram in Figure [5] shows the 
contribution of the background channels normalized to 
the known branching ratios [ifil Il7l | . It can be seen that 
the agreement between the background estimation and 
data is good near the K^K^ peak, indicating that the 
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FIG. 2: The K% momentum distribution. The dots with error 
bars are data, the dark shaded histogram is from K% mass 
side band events, and the light shaded histogram is from the 
Monte Carlo simulated backgrounds. The curve shown in the 
plot is the best fit of the data. 



estimation of the background under the KgK^ peak is 
reliable. 

Under SU(3) symmetry, KgK^ production via vir- 
tual photon annihilation is forbidden. This is checked 
by applying the same selection criteria to the data sam- 
ple taken below the tp(2S) peak, at -^i = 3.650 GeV, 
with an integrated luminosity of about one-third of that 
at the ip{2S) peak. Figure O shows the Kg momentum 
spectrum of the selected events; the events in the sig- 
nal region agree well with expectation from the Kg mass 
side band events. Taking all four candidates with mo- 
mentum between 1.7 and 1.9 GeV /c as signal, the upper 
limit of the production cross section at the 90% C. L. is 
measured to be a < 5.9 pb. The background from the 
continuum contribution is thus neglected in the folowing 
analysis since no evidence for KgK^ production via the 
virtual photon process is observed. 

The Kg momentum spectrum of the selected events 
is fitted from 1.45 to 2.0 GeV/c with a Gaussian distri- 
bution for the signal and an exponential for the back- 
ground using the unbinned maximum likelihood method. 
The result is shown in Figure |21 The backgrounds from 
the Kg mass side bands and the simulated background 
channels are also shown, and they agree with the fitted 
background reasonably well near the signal region, con- 
sidering the uncertainties in the global normalization of 
the background channels. The peak Kg momentum is 
(1775.0 ± 3.3) MeV/c, which is in good agreement with 
the Monte Carlo expectation. The momentum resolu- 
tion also agrees with the Monte Carlo simulation well. 
The fit yields n°^* = 156 ± 14 events, and the efficiency 
for detecting ■iP{2S) K^Kl, with K^ tt+tt" is 
epjc = (41.59 ± 0.48)% from the Monte Carlo Simula- 
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Kg Momentum (GeV/c) 

FIG. 3: The Kg momentum distribution for data at ^/s — 
3.65 GeV. The dots with error bars are data; the shaded his- 
togram is for Kg mass side band events. The curve shown in 
the plot is the best fit of the data, while the dotted histogram 
is the expected shape of a KgKl signal as determined by 
Monte Carlo simulation (not normalized). 



tion, where the error is due to the limited statistics of 
the Monte Carlo sample. 

Due to the long decay length of the high momentum 
Kg particles and the trigger requirement for hits in the 
VC, the trigger efficiency of KgKl events is very differ- 
ent than for normal hadronic events. Since the trigger 
system is not included in the Monte Carlo simulation, 
the trigger efficiency is measured using K^K^ events by 
comparing the number of events beyond and within the 
outer radius of the VC with what would be expected 
for an exponential decay, which yields a trigger efficiency 
of etrig — (76.0 ± 1.8)%. The systematic error in the 
branching ratio measurement from this source and all 
other sources are listed in Table. ^ 

The efficiency of the secondary vertex finding, e2nd, is 
studied using J/tp ^ K*{892)K + c.c. events. It is found 
that the Monte Carlo simulates data fairly well. Ex- 
trapolating the difference between data and Monte Carlo 
simulation to the Kg momentum range under study and 
correcting by the polar angle dependence of the efficiency, 
a correction factor of (98.1 ± 4.0)% is obtained to the 
Monte Carlo efficiency. 

The effect of the cut on the total energy of the photon 
candidates is checked with J/-)/; — > KgK^ events. The 
efficiency difference between data and Monte Carlo sim- 
ulation is measured to be 0.99 ± 0.01; data is slightly 
lower than Monte Carlo simulation. No correction to the 
final efficiency is made, and 2% is taken as the systematic 
error on the efficiency associated with this cut. 

The simulation of the tracking efficiency agrees with 
data within 1-2% for each charged track as measured 



TABLE I: Summary of systematic errors. 



oourcc 


Systematic errors (%) 


MO statistics 


1.2 


Trigger efficiency 


2.4 


Secondary vertex finding 


4.1 


cut 


2 


MDC tracking 


4 


Fit range 


2.4 


Background shape 


3.0 




5 


B{K% TT+Tf-) 


0.4 


Total systematic error 


9.2 



using channels like J/-0 ^ AA and ip{2S) tt+tt" J/-0, 
J/ip — > fi'^ iJ,~ . The systematic error for the channel of 
interest is taken conservatively as 4% . 

The Monte Carlo simulated mass resolution and mo- 
mentum resolution of the Kg agree with those deter- 
mined from data within the statistical uncertainties. The 
cut on the Kg mass distributions at 2 standard deviations 
introduces a very small systematic bias and is neglected. 

Varying the lower and upper bounds of the fitting 
range results in a 2.4% change in the number of the 
events; using a second order polynomial for the back- 
ground parametrization causes a 3% change in the num- 
ber of events. These are taken as systematic errors. The 
systematic error in the total number of ip{2S) events, 
^ip{2S)i which is measured using inclusive hadrons in the 
same way as in Ref. IJ^, is taken as 5%. The systematic 
error on the branching ratio B{K'g — s- tt'^'tt") is obtained 
from the PDG ij] directly. 

Figure 0] shows the cosine of the K^ polar angle for 
KgKl events from ip{2S) decays; agreement between 
data and Monte Carlo simulation is observed. This 
distribution is also checked with a larger sample from 
J/tp — > KgK^, where the Monte Carlo simulation agrees 
with data very well. This indicates that the angular dis- 
tribution used in the Monte Carlo generator is correct. 

The branching ratio of ip{2S) KgK^ is calculated 
with 

B{tp(2b) KgK^) = — — — g — — . 

N-4,(2S)B(K^ ^ TT+TT ) 

Using numbers from above (listed in Table. ITl)l . one ob- 
tains 

B{Tpi2S) KlKl) = (5.24 ± 0.47 ± 0.48) x 10-^ 

where the first error is statistical and the second system- 
atic. 

Comparing with the corresponding branching ratio for 
J/i) K^Kl from BESII ((1.82 ± 0.14) x lO"'') [3, 
one gets 

_ BmS)^K%Kl) _ 

- B{ Jli,-.K%Kl) - (2^-^ ^ 



5 




FIG. 4: Distribution of the cosine of the Kg polar angle of 
K%K% events from ip{2S) decays. Dots with error bars are 
data, and the histogram is the Monte Carlo simulation. 



TABLE II: Numbers used in the branching ratio calculation 
and the branching ratio result. 



quantity 


Value 




156 ± 14 


EMC (%) 


41.59 ±0.48 


£trig (%) 


76.0 ± 1.8 


£2nd (%) 


98.1 ±4.0 


iV^(2S)(10') 


14.0 ±0.7 


B{K°s -K+TT-) 


0.6860 ± 0.0027 fl7] 




5.24 ±0.47 ±0.48 



where the common errors in the J /ijj and i^{2S) analyses 
are removed in the calculation. The result indicates that 
ip{2S) decays are enhanced by more than 4cr relative to 
the "12% rule" expected from pQCD, while for almost 
all other channels where the deviations from the "12% 
rule" are observed, ip{2S) decays are suppressed. 

The branching ratio of KgR^, together with branch- 
ing ratios of '0(25') — > tt+tt" and '0(25') — > K^K~, can 
be used to extract the relative phase between the three- 
gluon and the one-photon annihilation amplitudes of the 
ipi^S) decays to pseudoscalar meson pairs. It is found 
that a relative phase around ±90° can explain the exper- 
imental results [111] . 

In summary, the flavor SU(3) breaking process KgKl 
is observed for the first time in "0(25) decays with the 
BESII i^{2S) data sample, and the branching ratio is de- 



termined to be i3(0(2S') ^ K^Kl) = (5.24 ± 0.47 ± 
0.48) X 10~^. Compared with the branching ratio of 
J/ip KgK^, 0(25) decays are enhanced relative to the 
"12%" pQCD prediction. The phases of the three-gluon 
and the one-photon annihilation amplitudes of 7/1(25) de- 
cays to pseudoscalar meson pairs are found to be nearly 
orthogonal. 
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